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Abstract 

Correlation functions, constructed from directional projections of the relative ve- 
locities of fragments, are used to determine the shape of the breakup volume in 
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coordinate space. For central collisions of Xe + Sn at 50 MeV per nucleon in- 
cident energy, measured with the 47r multi-detector INDRA at GSI, a prolate shape 
aligned along the beam direction with an axis ratio of 1:0.7 is deduced. The sensitiv- 
ity of the method is discussed in comparison with conventional fragment-fragment 
velocity correlations. 

Key words: Heavy ion collisions, multifragmentation, breakup state, correlation 
functions 
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A recent study of central ^^^Xe+'^'^^Sn and ^^''Au+^^'^Au collisions at bom- 
barding energies between 50 and 100 MeV per nucleon has shown that a good 
statistical description of the measured fragment yields and kinetic energies can 
be obtained provided that a prolate source deformation and a superimposed 
collective motion are included pj. The experimental data had been collected 
with the 47r INDRA multidetector [2] at the GSI laboratory. The statistical 
model employed in this study was the Metropolis Multifragmentation Monte- 
Carlo (MMMC) model [3] which had been extended to non-spherical (NS) 
sources P], a version referred to in the following as MMMC-NS model. The 
MMMC Statistical Model is based on the microcanonical ensemble and has 
found many applications in nuclear multifragmentation (see, e.g., [311511617118] ). 

A prolate deformation of the emitting source along the direction of the incident 
beam was indicated by the observed anisotropies of the fragment production. 
The element spectra were found to extend to larger atomic numbers Z at 
forward and backward emission angles than at sideward angles. The largest 
fragment within an event is preferentially emitted toward forward or back- 
ward angles. If emitted sidewards, its mean Z decreases, e.g., from Z ^ 18 
at 6'cm = 0° to Z ^ 14 at 6'cm = 90° for the 1% most central colhsions of 
^^^Xe+'^'^^Sn at 50 MeV per nucleon, selected on the basis of the measured 
charged-particle multiplicity or the transverse energy of light charged parti- 
cles pp. 

In the model description, an important role is played by the Coulomb inter- 
action which favours large separations between heavy fragments in order to 
minimize the Coulomb energy. Heavy fragments are, therefore, preferentially 
placed in the tips of a prolate source. Coulomb repulsion and the superimposed 
radial flow transform these spatial correlations into correlations in momentum 
space which produces the observed maxima in the yields and kinetic energies 
of the heaviest fragments at forward and backward directions. The orientation 
along the beam axis is clearly of a dynamical origin. The question, therefore, 
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arises whether the deformation of the model source, besides the superimposed 
flow, is mainly required for simulating the observed anisotropics in momentum 
space or whether it actually reflects the source shape at breakup as caused by 
the reaction dynamics. 

It is the aim of the present work to confirm the indicated source deformations 
in coordinate space by applying interferometric methods to the same data, 
here at first only for the case of ^^^Xe+'^'^^Sn at 50 MeV per nucleon. Interfer- 
ometry has become a standard tool for investigating the space-time properties 
of the breakup state in heavy- ion reactions [9lll0fll] . The spatial dimensions 
or space and time are usually separated by exploiting the effects of quantum 
statistics in proton-proton, neutron-neutron, or pion-pion correlation functions 
|10|12|13|14|15p6|| . Particle pairs are selected according to the orientation of 
their relative velocity with respect to the chosen coordinate axes, and separate 
correlation functions are generated as, e.g., for longitudinal and transverse ori- 
entations. In fragmentation reactions, the mutual Coulomb repulsion between 
fragment pairs has been used to derive time scales for the emission process 
from fragment-fragment correlation functions PT|18|19|2UII21] . Glasmacher et 
al. have shown that also information on the source geometry can be obtained 
if directional cuts are applied and correlation functions are generated for frag- 
ment pairs with relative velocities parallel or perpendicular to the symmetry 
axis [22123] • The dependence on the source geometry was small, even though 
nuclei with the unusual shapes of fiat disks or toroids were included in these 
studies. It will be shown in the following that the sensitivity to the source 
geometry is significantly enhanced with the proposed new kind of correlation 
functions of projections. 

The new correlation functions are constructed from the relative velocities of 
fragment pairs by not using their modulus, as for conventional fragment- 
fragment correlations, but rather their longitudinal and transverse compo- 
nents with respect to the direction of orientation which here is the beam axis. 
To correct, in first order, for the variation of the Coulomb repulsion with the 
fragment Z, the reduced velocity fred = "^^rei/ \J Zi + Zj is used, where frei, Zi 
and Zj are the relative velocity and the atomic numbers of the two fragments, 
respectively [21]. With the help of model calculations, it will be shown that 
these projected correlations are particularly sensitive to a deformation of the 
source. This will be done in comparison to conventional correlation functions 
generated not with directional cuts but with directional weights, proportional 
to sin^6' and cos^6' where 6 is the polar angle of the relative velocity vector 
with respect to the beam axis [16]. The method will then be applied to the 
experimental data for the ^^^Xe-|-'^'^*Sn reaction at 50 MeV per nucleon. 

Experimental details of these measurements, performed at the GSI laboratory 
in 1998 and 1999, and of the analysis and calibration procedures may be 
found in [T]|25|26ll27j . Natural Sn targets of areal density 1.05 mg/cm^ were 
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Fig. 1. Fragment-fragment velocity correlations as a function of the reduced velocity 
fred (in units of W~^c) of fragment pairs with 3 < Z <30 (top and middle panels) 
and of the two fragments with the largest Z within an event (bottom panel) for 
central collisions of ^^^Xe+^'^^Sn at 50 MeV per nucleon (dots). Solid, dashed and 
dotted lines represent the filtered predictions of the MMMC-NS model for the pro- 
late, oblate and spherical sources, respectively. The uncorrelated denominator was 
generated after azimuthal event alignment for the top panel and after polar event 
alignment for the middle and bottom panels (see text). The interval 50 — 100- 10~^c 
is chosen for normalization. 



bombarded with ^^^Xe beams delivered by the SIS heavy-ion synchrotron. 
Central impact parameters were selected by requiring that the total transverse 
energy E±i2 of light charged particles {Z = 1,2) exceeded 637 MeV [1]. This 
corresponds to an impact-parameter range 6/&max < 0.1, according to the 
geometrical prescription of [28] . 

Standard fragment-fragment velocity correlations constructed for this event 
class are shown in Fig. [H They were obtained with the technique of event 
mixing, and the interval 50 — 100 • lO^^c was chosen for normalization (c is the 
velocity of light). In the top panel, the uncorrelated denominator was gener- 
ated after rotating all events, azimuthally with respect to the beam axis, into 
a common reaction plane in order to suppress artificial flow effects |29]. The 
orientation of an event was determined from that of the longest eigenvector 
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of the kinetic-energy tensor, built from all detected charged particles. A flat 
correlation function at large f red is obtained by aligning all events with respect 
to each other before generating the uncorrelated denominator, e.g. by rotating 
their longest eigenvectors into the beam direction (Fig.[T], middle panel). With 
this polar event alignment, the distribution of mixed pairs extends to slightly 
larger v^ed, leading to a reduction of 1 + R with respect to the case of azimuthal 
alignment only (top panel). The simpler method of azimuthal event alignment 
has, nevertheless, been used in the analysis described below (Figs. [2lll]). 

All correlation functions are characterized by a strong suppression of small rel- 
ative velocities, with 1 + R falling below 0.5 for v^cd < 14-10~^c. Similarly large 
suppressions have been observed for a variety of heavy-ion and light-ion in- 
duced reactions and were consistently interpreted as indicating short breakup 
times of the order of 50 - 100 fm/c PUII^(.3 0.31J. The correlation function of 
the pairs with largest Z in each event (Fig. [H bottom panel) exhibits a con- 
siderable enhancement around 30 ■ 10~^c. With mean values (Zi) = 13 for the 
largest and (Z2) = 8 for the second largest fragment [1], the bump corresponds 
to relative velocities of ~ 4 cm/ns, a value typically observed for large-angle 
fragment-fragment correlations 032] • The distribution of relative angles be- 
tween the two largest fragments in the center-of-mass system is indeed wide 
and has a mean value of ~ 100°. 

The MMMC-NS model [4j has been used as an event generator for constructing 
correlation functions to be compared with the experimental data. For the 
global source parameters the previously obtained values of total charge Z = 79 
and mass A = 188, excitation energy = 6 MeV per nucleon, and collective 
flow 2 MeV per nucleon with a profile parameter ttcou = 2 were used [1] . Only 
the shape parameters were varied, and three different source elongations were 
chosen, a prolate source with axis ratio 1:0.70, an oblate source with axis 
ratio 1:1.67, and a spherical source. The produced MMMC-NS events have 
been filtered using a software replica of the experimental apparatus that takes 
into account the main properties of the INDRA multidetector as, e.g., the 
exact geometry of the detectors, the energy thresholds and limits for charged- 
particle detection and identification, and the effects of the multihit treatment 
in the off-line analysis. 

The standard correlation functions depend very weakly on the source shape 
(Fig. [1]). Only the correlation function for the two largest fragments exhibits 
slightly different maxima for deformed sources whose normalization, however, 
in the absence of an uncorrelated plateau region, is less certain (Fig. [U bot- 
tom panel). The depression at small fred is well reproduced, even though 
the predicted rise around fred = 15 ■ 10~^c is steeper than found experimen- 
tally. Smoother correlation functions could probably be generated and a better 
agreement reached with finite emission times and by allowing the source pa- 
rameters to fluctuate (cf. Ref. [33]). 
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Fig. 2. Directional correlation functions obtained from filtered predictions of the 
MMMC-NS model for fragment pairs with 3 < Z < 30 (left panels) and for the 
two fragments with largest Z within each event (right panels). The top and bottom 
panels show the transverse and longitudinal correlation functions, respectively, ob- 
tained with the method of directional weights [16j. Solid, dashed and dotted lines 
represent the predictions for the prolate, oblate and spherical sources, respectively. 
The interval 30 — 100 • lO^^c is chosen for normalization. 

The dependence of conventional directional fragment-fragment correlation func- 
tions on the shape of the emitting source in coordinate space is illustrated in 
Fig. [21 Longitudinal and transverse correlation functions were generated with 
the method of directional weights [T6], i.e. by accumulating the numerators 
and denominators with weight functions cos^^ and sin^6', respectively. Their 
dependence on the source shape is visible but weak, similar to the results 
reported in [22|23] . 

A much stronger sensitivity to the deformation is exhibited by the correlation 
functions generated separately for the longitudinal and transverse projections 
of the relative velocity fred (Fig. [3]). Note that for the projections shown here 
and below different intervals were chosen for normalization (see captions). 
Since individual components can become small even if the modulus of fred 
is large, these correlation functions do not approach zero at small values of 
the projected Wred- Their magnitude there is sensitive to correlations of the 
orientations of the relative- velocity vectors which, apparently, depend strongly 
on the source shape. The correlations of the transverse projections at small 
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Fig. 3. Correlation functions constructed from the transverse (left panels) and lon- 
gitudinal (right panels) projections of the reduced relative velocity of fragment pairs 
with 3 < Z < 30 (top panels) and of the two fragments with largest Z (bottom 
panels) for central ^^^Xe+^'^'Sn collisions at 50 MeV per nucleon (symbols). Solid, 
dashed and dotted lines represent the filtered predictions of the MMMC-NS model 
for the prolate, oblate and spherical sources, respectively. The dashed-dotted lines 
(closely following the dotted lines at small Ured) represent the spherical source with 
ellipsoidal flow (see Ref. [1]). The normalization was performed, individually for 
each projection, within 30 — 40 • 10~^c in the top and 20 — 30 • 10~'^c in the bottom 
panels, respectively. 



■^redX; paitlcular, change significantly with the source geometry. As will 
be shown and discussed below, the sensitivity to radial flow is large as well. 
With the flow value 2 MeV per nucleon determined from the fragment kinetic 
energies, a good agreement with the experimental data for ^^^Xe+'^'^^Sn is 
reached by assuming a prolate deformation. This is valid for the correlation 
functions generated from all pairs with j > 2 but, in particular, also for the 
correlation functions for the two largest fragments of each partition (Fig. [3l 
bottom) . The comparison thus confirms the conclusions drawn from the study 
of the fragment anisotropics in yields and kinetic energies [1]. 

Projected correlation functions generated with the MMMC model for central 
^^^Xe+'^'^^Sn collisions with the same input parameters but with the collective 
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Fig. 4. Correlation functions constructed from the transverse (left panel) and longi- 
tudinal (right panel) projections of the reduced relative velocity of fragment pairs 
with 3 < Z < 30 for central i29xe+"^*Sn colhsions at 50 MeV per nucleon (sym- 
bols, same data as in Fig. O top panels). Solid, dashed and dotted lines repre- 
sent the predictions of the MMMC-NS model without radial flow for the prolate, 
oblate and spherical sources, respectively. The normalizations were performed within 
30 - 40 • IQ-^c. 



radial flow set to zero are shown in Fig. HI Without flow, their rise at small Wred 
is much steeper, notably for the transverse projections, and their differences 
for the three considered source shapes are much smaller. It is evident that the 
experimental data for ^^^Xe+'^'^^Sn are qualitatively different from this set of 
correlation functions and cannot be reproduced without assuming collective 
contributions to the particle and fragment energies. 

The deduced sensitivities to the source shape in coordinate space were inde- 
pendently verified by reproducing the experimental correlation functions with 
a simpler source model [31]. Fragments were sampled according to the exper- 
imental yield distribution and placed randomly without overlap into a pre- 
determined ellipsoidal volume with the transverse {Rx,Ry) and longitudinal 
[Rz) source radii being varied between 6 fm and 30 fm. The parameters de- 
scribing the thermal and anisotropic collective motions of the fragments were 
chosen, individually as a function of Z, so as to best reproduce the measured 
kinetic energies. The temporal evolution after freeze-out was modeled with 
N-body Coulomb trajectory calculations. The asymptotic velocities were used 
to construct both, weighted and projected, correlation functions. The result 
of a comparison with the experimental data, for intermediate-mass fragments 
with 5 < Z < 7, is shown in Fig. [5] in the form of the distributions in the 
plane of the transverse and longitudinal extensions, representing the degree of 
agreement obtained. 

The different sensitivities of the two types of correlation functions are imme- 
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Fig. 5. (color online) Contour plots (in units of the standard deviation a) of the 
two-dimensional distributions in the plane of the transverse (xy) and longitudinal 
(z) extensions of the source, as obtained from the comparison of the experimental 
weighted (left panel) and projected (right panel) correlation functions, for fragments 
with Z = 5 — 7, with the predictions of the classical deformed-source model for 
^^SXe+'^^^Sn at 50 MeV per nucleon. 

diately apparent. With directional weights, the minima of the distribution 
extend over a broad band of radius values, all corresponding to approximately 
the same volume but to different shapes. For the correlation functions con- 
structed from the projections (right panel of Fig. [5]), the minima are concen- 
trated in a region of prolate source shapes with axis ratios of 1 : 0.6 ± 0.1. 
Other observations were also found to be similar to those made with the 
more complex MMMC model. In particular, the requirement of a finite radial 
flow for obtaining the shape dependence of the projected correlation functions 
was confirmed. Radial flow leads to a correlation between the locations of a 
particle in coordinate space and in momentum space which makes its role 
here understandable. However, flow by itself is not sufficient (cf. dotted and 
dashed-dotted lines in Fig. [3]). The projected correlation functions become 
significantly different only in the presence of deformations in coordinate space 
(full and dashed-dotted lines in Fig. [3]). 

The volumes obtained with this method are large. A mean radius of 15 fm 
(Fig. El left panel) corresponds to about twice the radius of a system of A = 188 
nucleons at normal density, slightly larger than what is usually assumed in 
applications of the MMMC model (6 times the normal nuclear volume in 
Refs. [1]) or what is obtained with other methods [35|l36] . A very similar vol- 
ume results from the transverse and longitudinal radii of 13 fm and 22 fm, 
respectively, favoured by the method of projections (Fig. right panel). Con- 
sidering the finite range of Z used for this comparison, this is perhaps not 
unreasonable. It is, furthermore, known that large variations of the volume 
are necessary in order to change fragment-fragment correlation functions sig- 
nificantly [37]. An additional contribution to the apparent source size may 
also come from the finite time scale of the emission process even though it is 
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expected to be rather short, perhaps 50 - 100 fm/c (see, e.g., Ref. [23f38] ). Es- 
timates for the increase of the mean source radius, obtained from simulations 
with the above-mentioned simpler source model using finite emission times, 
are of the order of 1 - 2 fm. They are non-negligible regarding the correspond- 
ing breakup volume but they are unlikely to change the parameters describing 
its anisotropy. 

In summary, a new method based on correlation functions of relative-velocity 
projections has been presented and shown to be sensitive to source defor- 
mations in the presence of flow. It has been applied to the data for central 
collisions of ^^^Xe+'^'^^Sn at 50 MeV per nucleon, and an expanded source with 
a prolate elongation in the direction of the incident beam has been deduced. 
The axis ratio 1 : 0.7 used for the quantitative comparison with the MMMC- 
NS model predictions has led to very satisfactory results, thus supporting the 
assumption of a non-spherical breakup source with a deformation of this mag- 
nitude for the statistical description of the fragmentation channels with this 
model [Tj. A similar deformation of 1 : 0.6 ±0.1 has resulted from the analysis 
with a simpler source model, restricted to the Z = 5-7 range of fragments. 

The dynamical origin of this deformation was not a subject of the present study 
but is most likely found in the incomplete mutual stopping of the incident ions, 
even in the most violent collisions associated with the largest transverse en- 
ergies or particle multiplicities. Stopping is incomplete even in the heavier 
^^''Au+^^^Au system for which an excitation function of stopping with a max- 
imum at several hundred MeV per nucleon has been established [H^PP] . At the 
present energy of 50 MeV per nucleon and below, the observed mass hierarchy 
of the longitudinal fragment velocities provides further evidence for dynamical 
correlations between the entrance and exit channels [H]. In coordinate space, 
the larger longitudinal momenta will cause the source shape to be elongated 
at the time when the breakup density has been reached. Larger fragments, 
presumably carrying most of the remaining longitudinal momentum, will be 
near the tips of the elongated source, an effect which in the static MMMC- 
NS description is generated by the Coulomb force alone. Coulomb effects are 
probably less dominant in the dynamical situation which explains why the 
deformed-source model, by its construction, permits exploring the statistical 
nature of multi-fragment breakups in the presence of dynamical constraints 
which are not part of the model. Source deformations and anisotropic emis- 
sions are generic features of central collisions at intermediate energies, a fact 
already quite well established but supported in a new way with the present 
method. Since the radial flow is increasing it should be possible to extend 
these studies to higher bombarding energies. 

This work was supported by the European Community under contract No. 
ERBFMGECT950083. 
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